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Abstract. Spectral enhancement - which aims to undo spectral broadening 
— leads to integral equations which are ill-posed and require special regular- 
isation techniques for their solution. Even when an optimal regularisation 
technique is used, however, the errors in the solution - which originate in data 
approximation errors — can be substantial and it is important to have good 
bounds for these errors in order to select appropriate enhancement methods. A 
discussion of the causes and nature of broadening provides regularity or source 
conditions which are required to obtain bounds for the regularised solution of 
the spectral enhancement problem. The source conditions do only in special 
cases satisfy the requirements of the standard convergence theory for ill-posed 
problems. Instead we have to use variable Hilbert scales and their interpola- 
tion inequalities to get error bounds. The error bounds in this case turn out 
to be of the form 0(e 1— 'Jf 6 )) where e is the data error and 77(e) is a function 
which tends to zero when e tends to zero. The approach is demonstrated with 
the Eddington correction formula and applied to a new spectral reconstruction 
technique for Voigt spectra. In this case 17(e) = 0(l/V]loge|) is found. 

1. Introduction 

One of the computational challenges in spectroscopy is the separation of overlap- 
ping spectral lines. This separation can be achieved by computationally narrowing 
the spectral lines and thus enhancing the resolution or correcting the spectrum. 
The class of methods of resolution enhancement considered here is based on the 
solution of linear Frcdholm integral equations of the first kind using observed data 
for the right hand side. The basic approach was first analysed in [3] but it goes back 
in principle to work by Stokes [UJ . The effect of data errors has to be analysed 
carefully, especially since the enhancement problem is ill-posed. This analysis is 
performed in the following using variable Hilbert scales [531 HI] ■ A more traditional 
error analysis which can be found in |20j is not directly applicable here as the source 
conditions are non-standard. However, in contrast to many other ill-posed prob- 
lems, here the underlying physical model does suggest specific source conditions. If 
/ is the enhanced spectrum and f a an (optimal order) regularised approximation 
of / then bounds of the form 

\\f - f a \\ < e 1 ^ 

are found where e is the residual of f a . In the classical case the 77(e) is constant, in 
contrast it is shown here that this exponent slowly decreases to zero with e — > 0. 

A new enhancement method based on Lorentz kernels for Voigt spectra is shown 
to provide good performance compared to more traditional methods like the Ed- 
dington correction as it capitalises more on the smoothness of the data and does not 
require any advanced knowledge of the proportions of the Gaussian and Lorentzian 
components in the Voigt spectrum. If a spectrum contains a Gaussian component 

the error bound is of order 0(e 1 ~ c '^ log£| ) and the convergence rate thus grows 
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with £->0. For very small e one can find very close to O(e) convergence, however, 
this depends on the level of enhancement required. Experiments show that this 
method leads to a reduction of linewidth of more than a factor of two in the case 
of a 5% data error. 

In the remaining parts of this section a brief review of broadening mechanisms 
are given, in addition to a short discussion of a least squares method to determine 
the location and strength of spectral lines. In section 2 we present the integral equa- 
tion framework for resolution enhancement and illustrate this with the Eddington 
correction formula and Stokes correction by partial Gaussian deconvolution. In 
section 3 the method using Lorentz deconvolution for Gaussian and Voigt spectra 
is discussed in terms of the errors. Section 4 then provides some demonstrations 
of the enhancement properties of this Lorentz deconvolution which in particular 
illustrates the broadening effects of noise and regularisation. In the concluding 
section 5 related and open problems are considered. 

1.1. Models of spectra and broadening. In the natural sciences, a spectrum is 
a distribution of photon counts over energy or frequency. Since Fraunhofer's work 
in 1814 it is well known that this distribution is concentrated along lines, both for 
emission and absorption spectra. The existence of these spectral lines was later 
confirmed by quantum mechanics. Their importance is due to the fact that they 
provide information about the energy levels of the electrons and thus insights into 
the structure and composition of the originating substrate. Spectroscopy has been 
for a long time one of the most important tools in experimental science. A simple 
model for a spectrum based on the Fraunhofer spectral lines would consist of a 
probability measure with discrete support. 

Almost simultaneously with Fraunhofer's discovery it was realised that spectral 
lines have a non-zero width. This broadening originates from many different phys- 
ical effects and a discussion of spectral broadening can be found in a variety of 
different books and journals, see for example [HI El 1121 EDI El SH El C31 HH1 EZJ- In 
order to get a basic idea we review some of the most important mechanisms here. 

A first type of broadening, termed natural broadening, occurs because the time 
of the transition between the two energy levels is finite. The spectral lines which 
have only been broadened by this type have a Lorentzian shape, i.e., have peaks of 
the form 1/(1 + x 2 / s 2 ) where s is a width parameter. Usually natural broadening 
leads to very narrow lines. Much larger than natural broadening is usually Doppler 
broadening which occurs because the emitting (or absorbing) particles are in con- 
stant thermal motion which leads to a Doppler effect which shifts the energies of 
the photons. The shape of spectral lines which only have been Doppler broadened 
are Gaussian. While the width of the Doppler broadened lines is proportional to 
the energy we will neglect this here and assume a constant width approximation. 
Neighbouring particles to the electrons emitting or absorbing the photons produce 
a third kind of broadening, the pressure broadening. One can show that in the case 
where only pressure broadening occurs the spectral lines are Lorentzian. Further 
broadening originates in the instrumentation and even discretisation (or binning) 
of the spectrum produces a certain amount of broadening |14j . Finally the medium 
which the photons need to traverse before getting to the observer also produces 
some broadening. There are other effects which contribute to broadening and there 
are other distortions of spectra than broadening occurring. This includes spectral 
shifts and the occurrence of extra peaks, so-called satellites [LI] . 

A fairly general but simple broadening model would represent observed spectra 
as the effect of an integral operator on an underlying spectrum which might have 
been modified in other ways. Here this underlying spectrum u is assumed to be in 
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Li2(M) and so an observed spectrum g is of the form 

g(x) = / a(x,y)u(y)dy 



R 



with some kernel a which in the simplest case is assumed to be a convolution 
kernel, i.e. a(x, y) — ot{x — y) for some L2 function a. More generally, an observed 
spectrum is modeled as the image of a product of several broadening operators 
A\,..., A n , i.e., as g = A\ ■ ■ ■ A n u. In some cases, such a product can lead to a 
normal distribution because of the central limit theorem. Here we assume mostly 
that all the operators are convolutions and have Lorentzian or Gaussian shape (but 
different widths) . As the operators commute and the convolution of Lorentzians is 
a Lorentzian and of Gaussians is a Gaussian, respectively, it is found that a good 
model is given by the Voigt shape which consists of a convolution of a Gaussian with 
a Lorentzian. In the following we call the integral equation Au = g representing 
any kind of (linear) broadening the broadening equation. 

1.2. Fitting the lines. While immediately appealing, the inversion of the broad- 
ening equation Au = g is not feasible as it is typically severely ill-posed, the g has a 
substantial amount of observational error and u is typically not very smooth so that 
even a regularised solution cannot be expected to be a good approximation. Any 
feasible approximation will make use of the (approximate) Fraunhofcr line structure 
of the u. The simplest model assumes that u is a measure with discrete support 
and intensities Ui so that the broadening equation takes the form 



The determination of the Xi and U{ from some data gs with \\g — gg\\ < 5 can be 
done by minimising the least-squares objective function 



When the locations Xi of the spectral lines are known this amounts to a linear least 
squares problem. The determination of these locations, however, is a nonlinear 
problem. An interesting discussion of this problem from the perspective of Bayesian 
statistics can be found in [TO] . 

In |19j Golub and Pereya discuss the variable projection method for the solution 
of the nonlinear problem above in the case of a finite number of non-zero itj . Rather 
than minimising the squared residual they first solve for the linear parameters itj 
explicitly such that u = A + {x)gs where u = (ui, . . . ,u n ). They then use a non- 
linear (typically Gauss-Newton) method to solve for the locations x — (xi, . . . ,x n ) 
by minimising the functional || A(x)A + (x)gs — g$\\. In a recent paper |39j the au- 
thors discuss the application of this method to spectroscopic problems and consider 
reasons for the success of the approach. They observe in particular superior numer- 
ical conditioning and convergence of the Gauss-Newton method compared to the 
original optimisation problem. 

An important condition required by the variable projection method is that the 
matrix A(x) has to have a fixed rank for x in some neighbourhood of the minimum 
of the variable projection functional. This condition may be difficult to fulfil when 
one has two components of x which are very close. As two coinciding Xj will reduce 
the rank of A(x), the neighbourhood where the rank condition holds can be very 
small. It would certainly be difficult to find initial conditions for the Gauss-Newton 
iteration which are in a neighbourhood of the exact solution. 
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When the spectral lines are well separated then the variable projection method 
works very well. This is for example the case where the baseline condition in which 
case the functions a(-,Xi) have non-overlapping supports (at least numerically). 
It follows that the a(-,Xj) are pair- wise orthogonal, good starting values can be 
obtained and the rank condition can be maintained. A similarly favourable situation 
occurs if the Rayleigh condition holds. This motivates the development of methods 
which are able to enhance the spectrum so that the enhanced spectral lines are 
better separated. A discussion of these aspects from a statistical perspective can 
be found in pQ. 

2. Resolution enhancement 

2.1. The enhancement equation. The resolution enhancement procedures con- 
sidered here consist of algorithms which determine the enhanced spectrum f as 
a solution of an integral equation Bf — g from the observed spectrum gg which 
satisfies \\gs — g\\ < 5. The integral operator B is of the form 

(1) Bf(x)= f b(x,y)f(y)dy. 
The integral equation 

(2) Bf = g 

will be called the enhancement equation. The operator B is chosen such that the 
enhanced spectrum / has narrower lines than the original spectrum g. The main 
constraint in choosing B is that the enhancement equation should be solvable which 
means that g has to be in the range of B: 

(3) g e range (B). 

In the case where B is a convolution operator, the resolution enhancement is the 
Stokes correction formula |43j . The integral equation Ansatz for enhancement was 
introduced Allen, Gladney and Glarum in their ground-breaking paper [2 . A simple 
precursor to this type of enhancement is the Eddington correction formula |16l 1171 [8] 
for the enhancement of spectra with Gaussian peaks using differentiation. 

The careful choice of the operator B is essential to successful enhancement. Even 
if the range condition Q holds, the solution of the enhancement equation (J2]) may 
show poor resolution and contain a large error. This is due to the ill-posedness of 
the enhancement equation. It's solution will require some form of regularisation. 
When selecting B one has to trade-off the amount of enhancement achievable by 
B against the regularisation required for the solution of the enhancement equation. 
While the theory of resolution enhancement is based on the general theory for the 
solution of integral equations, there is one important difference: When solving inte- 
gral equations, the operator is given while for resolution enhancement, the operator 
B is chosen. In both cases, on needs to choose the regularisation method. 

There is a large literature on regularisors, a concise and short reference is still 
the book by Groetsch [20! ■ In this book, convergence rates of regularisors are given, 
provided that a source condition of the form g E range((_Bi?*) s ) holds for some in- 
teger s > 1 and where B* denotes as usual the adjoint of the operator B. Here we 
will use a more general theory based on variable Hilbert scale inequalities [23l [24] . 
This framework has since been used in j32j [3H [3TJ [33] . In the analysis literature, the 
variable Hilbert scale interpolation is called interpolation with a function parame- 
teiQ see, for example [37J [T2J [3H]- In the analysis of partial differential equations, a 
related generalised Holder inequality has been applied in [B] . Source conditions are 
very important in the analysis of convergence of regularisation and some newer work 
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which includes the application to nonlinear problems can be found in [15 , 45 , 31 , 26 . 
The recovery of / = B~ 1 g from gs is the main topic of the book [21] by Groetsch. 
The specific case of singular convolutions are covered in a paper by Sushkov [44 . 

In the following let Hg C £2(1^) denote the Hilbert space with the norm \\g\\B = 
||-B -1 <7|| and let be the Hilbert space with the norm \\g\\^ — (g,ip((BB*)~ 1 )g) 
where tp is a function on (a subset of) (0, 00) which is continuous and monotonically 
increasing. The operator ^((BB*) -1 )) is defined using the spectral theorem as 
in [24]. In the following f a will always denote a regularised solution of Bf = g. 
One then has the following general convergence theorem. 

Theorem 1. Let B : L2(R) — > be an injective, continuous linear operator 

for which BB* — XI is injective for < A < c$ for some cq > 0. Furthermore, let ip 
be a non-negative function which is monotonically increasing for arguments larger 
than 1/cq. Finally, let ^ be a non-negative function such that ^(ijj(\)) > A and ^ 
is monotonically increasing and concave for all arguments A > 1 /cq . 
If fa £ Hip satisfies 

(4) \\Bf a \y<C, and 

(5) ||B/ a - ff ||=e 
then 

(6) \\f-f a \\<e^((C+\\gU)ye^ 
for all f and g = Bf € . 

Proof. The functions 'J, tp together with the functions 9(A) = 1 and 0(A) = A 
satisfy the conditions of Theorem 1 in [22] which is a direct consequence of the 
interpolation inequalities in [23]|24] which, with — \\r\\B and ||r||g = ||r|| here 
takes the form 

IH| s <||r|| y/nHl/WrW 2 ), for all re IV 
Now let r = Bf a — g. As g = Bf on has 

|M|B = ||B/a-S/||B = ||/a-/||. 

Furthermore, by the triangle inequality one gets 

\\rh = \\Bf a -gU<\\Bf a \\ i ,+ \\g\\^ 
and, as ^ is monotonically increasing, it follows that 

H\\r\\l/\\r\\ 2 )<m(C+\\g\M 2 /ei). 
Inserting this in the interpolation inequality gives the claimed bound. □ 

This result can be interpreted as a variant of the Lax equivalence theorem. 
The conditions on f a are the stability condition H-B/qH^ < C and the consistency 
condition \\Bf a — g\\ — e. If ip is unknown one may take a stronger norm for 
stabilisation in a discrepancy method similar to the one discussed in [23 . For 
consistency one wants to make sure that e is small. This is achieved indirectly by 
controlling the size of \\Bf a — gs\\ and observing that 

\\Bf a -g\\ < \\Bf a -gs\\ + \\g-gs\\ 

by the triangle inequality. In the following we call any (approximate) enhancement 
f a which satisfies both conditions Q and Q a spectrum which has been stably 
enhanced with B. 
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2.2. The Eddington correction formula. This early and still popular approach 
to the enhancement of Gaussian spectra uses derivatives and is of the form 

, 9 {2) , .9 (4) 

f = 9 ~^ + - 8 ' 

see [TSJ El [5] . It has been observed in [5] that correction formulas of this type may 



be viewed as solutions of integral equations of the form discussed in section 2.1 We 
can thus apply theorem [T] to obtain an error bound for the Eddington correction. 
See also [36] for a discussion of their application in practice. Other procedures to 
spectral enhancement based on differentiation are discussed from the point of view 
of numerical differentiation in [3] . 

The A:-th order Eddington correction / is defined as 

3=0 J ' 

where </ 2j ) denotes the derivative of order 2j of g. The Eddington correction 
formula now fits into the integral equation framework for resolution enhancement 
with enhancement equation Bf = g and the enhancement operator B has a kernel 




K x >y) = - Z^^-y cos((j(x - y))du. 



In particular, for k = 1 one has 

and for k = 2 the kernel is of the form 

b(x, y) = ie - alx - yl cos(p(\x -y\ + 9) 

for some a, j3, 7 and 6. 
In the following, let 

a G (x,y) = 4^e- (x - y)2/2 



/2ir 

and let a spectrum g which has been broadened by a G be called a Gaussian spec- 
trum. In this case one has 

(8) g{x) = I a G {x 1 y)u{y)dy 



for some u G L2(IR). The Eddington correction formula have been designed to 
reduce some of the broadening produced by a G . 

A motivation for this particular formula comes from the convolution theorem as 

g(u>) = a G {uj)u{uj) 

where a G (uj) — exp(— w 2 /2) and g and u are the Fourier transforms of g and u 
respectively. By the Taylor theorem one then gets formally 

3=0 ■'' 

Truncating this expansion and using the fact that multiplication with uj 2 in the 
Fourier domain corresponds to taking —d 2 /dx 2 in the original domain gives the 
formula. 

The following lemma provides the expressions and some properties for the ip and 
^ which will be used to establish the error bound of the correction formula. 
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Lemma 1. Let B be the enhancement operato^for the k-th order Eddington cor- 
rection formula. Furthermore, lettf.(r]) be the k-th order Taylor polynomial for the 
exponential function for k > and tf. — for k < 0. Then 

(1) Any Gaussian spectrum g is in H^, the Hilbert space with the scalar product 
(g,g)i> = {g^{{BB*)~ 1 )g) and where 

V>(A) = exp(2^ 1 (\/A)), A>1. 

(2) The inverse 

*(r ) )=ip- 1 (r,)=t k (log(r l )/2) 2 , r, > 1 

is concave. 

Proof. (1) The £?-norm is by Parseval's theorem 

\\g\\l = \\B- l gf= [ t k {w 2 /2) 2 \g{u)\ 2 du. 

JR 

As ■0(t/ c (w 2 /2) 2 ) = exp(w 2 ) by definition one gets 

\\g\U, = {g,ip{{BB*)- l )g) = J^M^)\g{^)\ 2 du>. 
which is equal to 1 1 it || 2 if g is a Gaussian spectrum with 
g(x) = I a G {x,y)u(y) dy. 



It follows that \\g\\ih is a norm on the set of Gaussian spectra which provides 
a Hilbert space structure for this space. 
(2) As dtk(()/d( — ifc_i(C) one has d^(Q/d( = t k -i(C) and consequently 

d 2 ^ _ i ( i (\og{o\ k tk _^ | 1 ziog^y- 1 \ 



□ 



d( 2 2( 2 U! V 2 J K ~ L (fc-1) 
which is non-positive and so ^(C) is concave for C, > 1 



We now get the main theorem which provides bounds on how well one can 
evaluate the Eddingtion correction. 

Proposition 1. Let f a be a stably enhanced spectrum using B the k-th order Ed- 
dington enhancement for Gaussian spectra and = exp(2t^" 1 (-\/A)). Then there 
exists a C > independent of e such that 

(9) 11/ <Ce|log( £ )| fc . 

Proof. By theorem [T] and lemma [I] one has for 1/e > C + ||<7||i/>: 

Wf-fa\\<eiJ*ttC+\\9h) 2 /e 2 ) 

<et fc (21og(C+|| 5 y-2log(e)) 
<ee2 k (log(C+\\gU)-log(e)) k 
< ee4 k (-log(e)) k 
<Ce|log(e)| fc . 

as ffe(A) < e\ k for A > 1 □ 



Some times the inverse B 1 is called enhancement operator 



8 



MARKUS HEGLAND 



A consequence of this lemma is that the ill-posedness of the problem is really an 
issue for very high derivatives only. However, it is necessary to use regularisation 
nonetheless as otherwise the data errors would remove any advantage of the resolu- 
tion enhancement and typically render the so "enhanced" spectrum useless. Allen 
et al. [2] provide similar correction formulas to the Eddington formula for Lorentz 
spectra and also provide other correction formulas determining the coefficients in 
different ways, see also [25 . The analysis of the accuracy of so enhanced spectra 
can be analysed in exactly the same way as the Eddington formula. 

In order to compare the above error bound for the Eddington correction formula 
with the ones which we will obtain for other enhancement methods, one could 
restate it as 

11/ - /all < Ce^ 

where the exponent is 

,log|lo^_ 
I log(e) | 

The formula is valid asymptotically and we assume that < e < 1/e. It can be 
seen that the smallest exponent is now obtained for e = e~ e as 

Vmin = 1 — k/e 

and consequently 

11/ -/all <C^- k l\ 

It follows that for k = 1 , 2 one gets an error bound which is similar to the one 
obtained for an enhancement obtained through sharpening, see 22J. One can also 
get similar bounds for larger k a necessary condition on the error in this case, 
however, is 

lo^JogM 
I log(e) | 

and while first and second order Eddington corrections (with second and fourth 
derivatives) should work well even in the case of larger errors, but for higher order 
derivative corrections one does require smaller data errors. 

2.3. Stokes enhancement with a Gaussian kernel. By using Fourier trans- 
forms, Stokes |43j was able to introduce more general spectral correction formulas 
which amount to general deconvolutions. An example of such a formula would use 
a Gaussian kernel of the form 

1 / (x - v) 2 

(10) b(x,y)= -^exp' 



2ttk V 2k 



One can see that a resolution enhancement using this kernel reduces the width of a 
unit Gaussian spectral line from equation ^ from one to \/l — k 2 . The enhanced 
spectrum is again a Gaussian with no other local maxima and no local minima. 
While such an approach can be generalised to other than Gaussian spectra (see [2"2~P 
it does require the knowledge of the spectrum. As Gaussian spectral lines are very 
smooth, using this type of enhancement for less smooth non-Gaussian spectra will 
lead to meaningless results as the range condition is not satisfied in such a case. 

For the Gaussian case, however, one has the following result about the error of 
a regularised enhancement f a : 

Proposition 2. Let g be a Gaussian spectrum which has been enhanced by an 



operator with kernel b given in equation (10). Then the stable approximation f a 
satisfies the error bound: 

(11) 11/ -/all <C^\ 
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Proof. Using Fourier transforms and the Parseval equality one derives = A 1 ^ . 
As k G (0, 1) the inverse ^(r/) = = r) K is concave and the bound then follows 

from theorem [T] □ 

Note that in this case the source condition is a of a classical form and thus the 
error bound may also be obtained using methods from [20] . 

As the spectral enhancement reduces the width by a factor yl — k 2 it follows 
for example that a reduction of the width by a factor two is obtained by solving an 
integral equation of the first kind with error 0(e 1 ^ 4 ) if a stable method is used and 
e is the data error. 



3. Enhancing Voigt spectra with unknown line shape 

While it is known that many spectra are of Voigt type, i.e., they contain a 
mixture of Gaussian and Lorentz broadening it is often unknown, how much of both 
types are current in any particular spectrum. We will now present an enhancement 
procedure which utilises a Lorentz kernel for the enhancement of a Voigt spectrum. 

The enhancement equation Bf — g providing the enhancement is an integral 
equation with a Lorentz kernel of the form 

b(^x,y^j — — . 

' K7T 1 + (x — y) 2 1 K 2 

Thus Bf is again a convolution and the Fourier transform is 

b(uj) — exp(— 

The width parameter k has to be chosen similar to the width parameter for the 
Gaussian sharpening discussed in section |2.3| or the order of the Eddington correc- 
tion formula of section 12.21 In this choice one considers the trade-off between the 
enhancement obtained through the narrower lines in the spectra and the error from 
the solution of the integral equation. 

Before discussing the general case of a Voigt spectrum we provide a bound for 
the error of the Stokes correction with Lorentz kernel of a Gaussian spectrum. 

Lemma 2. Let B be the enhancement operator for the Stokes correction formula 
with a Lorentz kernel with width n. Then a Gaussian spectrum is in the space 
(based on B) with ip(^) — ex P((l°g(A)/(2K)) 2 ). Furthermore, the inverse ^(r)) = 
l if)~ 1 (rf) = exp(2Ky / log(?7)) is concave if k < \/2 or if n > | + \J (k/2) 2 — 1/2. 

Proof. We use the Fourier transforms of the kernel of BB* which is exp(— 2k|cj|) 
and AA* which is exp(— uj 2 ) and the Parseval equality to get i/j. 
The second derivative of ^ is then 

d 2 ^ _ K (2 log t/- 2vTog^K+ 1) e 2 ^°^ K 

d7 f 2r, 2 (log 77) 3 

and it follows that ^> is concave if 2 log rj — 2 y/log r\ k + 1 > 0. The conditions then 
follow directly. □ 

It then remains to apply Theorem [T] to get the following error bound: 

Proposition 3. The error of a stably computed enhancement f a of a Gaussian 
spectrum using the Stokes correction formula with a Lorentz kernel is bounded by 

WU-fW^^-^V^i 

for < e < eo and some cq > 0. 
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Proof. By theorem [T| and lemma [2] one has for e > and some C which satisfy 
1 



> C 



+ \\9h > V«/2+yW2)" 2 -l/2 



the bounds 



\\f-f a \\<e^((C+\\g\Ur/e*) 

< eexp( K y / log((C+|| 5 || v ,) 2 A 2 )) 

< eexp(2«;yi log e| 

< £ l-2«/\/i loge|_ 

□ 

Note that here re is not the width of the enhanced spectrum but a parameter 
which controls how much enhancement is done. Thus a larger re corresponds to 
more enhancement and re — to no enhancement. For example, if one has e ~ 1CP 3 
and re — 0.7 one gets an error of approximately 0(e 1 ^ 2 ). 

As discussed in the introduction many spectra have undergone broadening both 
with Gaussian and with Lorentz kernels. The resulting class of spectra are the 
Voigt spectra. We assume here that we know that a given spectrum is in this class, 
however we do not assume that we know how much each of the two components 
have contributed to the broadening. This is why we suggest a Stokes correction 
with Lorentzian kernels. 

Specifically, let the Lorentz kernel be 

a>L(x,y) 



with Fourier transform 

a L (uj) = exp(-\/2|w|). 

The Voigt spectrum (with mixing parameter 9) is then defined by its Fourier trans- 
form 

a v (u>) = a G (uj) e a L (w) 1-9 , 

and the kernel is thus 

2tt J r 

A Voigt spectrum is then of the form 

g(x) = J a v (x,y)u(y) dy 
for some u e L2(M) and < 6 < 1. One then has 

Lemma 3. Let B be the enhancement operator for the Stokes correction formula 
with a Lorentz kernel with width k. Then a Voigt spectrum with parameter 9 is in the 
space (based onB) with ip(X) = cxp (o (log(A) /(2k)) 2 + 78(1-6*) log(A) /(2k) 
Furthermore, the inverse ^ = ip is defined by 



^(rj) = cxp 



2 k 



(V2(l - 9) 2 + fllogfa) - V2(l - eju 



and is concave if k < ^/29 3 ^ 2 or if r\ > 770 for some rjo > 0. 
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Proof. We use the Fourier transforms of the kernel of BB* which is exp(— 2k|cj|) 

and AA* which is exp(— 9u> 2 — \/8(l — 9)\uj\) and the Parseval equality to get ip. 

With Civ) = n/IosC^) + b , a = 2k6~ 3 / 2 and b = 2(1 - 9) 2 9 one then has for the 

second derivative of W 

r- d 2 ^ nP a C('?) 

exp(v^(l -8))— = -__(2C(r ? ) 2 - aCfa) + 1). 



dr\ 2 4rj 2 ((ri) 3 

One gets convexity for 'J if 2£(?7) 2 — a£(rf) + 1 > which happens if k < \/2 9 3 / 2 , 
or for 77 > 770 and large enough 770. □ 

Then an application of Theorem [T] provides again an error bound: 

Proposition 4. The error of a stably computed enhancement f a of a Voigt spec- 
trum with width parameter 9 using the Stokes correction formula with a Lorentz 
kernel is bounded by 

\\f a ~ /II < e 1 - 2K /%/ e l 1 °g e l+( 1 - e ) 2 
for < e < eo and some e > 0. 

Proof. By theorem [T] lemma [3] and the monotonicity of ^S> one has for e > satis- 
fying e(C+ Ugll^) < 1 the bounds 



11/ - f a \\<e^{(C+\\9h?/<?) 
Wf-f a \\<eVW^) 

< eexp Q (y2(l-0) 2 + 4#|loge| - ^2(1 - 9) 



= e l-»7(e) 



where 



77(e) = 



V2k ( v /(l-0) 2 + 20|log(e)|-(l 



V2k 



|loge| 

2 



< 



V(l-0) 2 + 20|loge| + (l-0) 
2k 



V^j log e| + (1-9) 2 

As < e < 1 an upper bound for 77(e) will lead to an upper bound for the error. □ 

4. Enhancing a Gaussian peak 

We provide some simple experiments which show how resolution enhancement 
modifies a single Gaussian peak. In Figure [T]a Lorentz correction formula is applied 
with different values of the parameter k. Comparing the widths at height 0.5 one 
sees that for k ranging from \[2 to 4 one gets reductions of the widths between 
a factor of 1/2 to almost 1/5. Note that resolution enhancement comes at a cost 
which grows with k in the sense that side bands start to occur. From the plot it 
appears that the peaks of the side bands are at the level of the original (unenhanced) 
spectrum but can be negative. In Figure [2] one sees how regularisation (using 
the source condition) does further distort the peak. In this case we choose k = 
2. One can clearly see the oscillations and the broadening which are caused by 
regularisation. Finally, Figure [3] considers the same regularisation methods (except 
for the case of the regularisation parameter a — which gives much larger errors) . 
Here a data error of 5% has been included. The effect of the error onto the enhanced 
signals is that mainly the oscillations away from the centre are strongly affected by 
data error, especially for the case of a small regularisation parameter. 
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Figure 1 . Lorentzian correction of a Gaussian 



5. Conclusion 

From the physics of spectral broadening one obtains the broadening equation 
Au = g. For various reasons including the severe ill-posedness of the equations, 
the fact that A might not be known, and that u might not be sufficiently smooth, 
the solution of Au = g is typically not feasible. However, this equation provides 
a regularity or source condition for the solution of enhancement equations Bf = 
g which are essential for obtaining error bounds or regularisation methods. As 
typically A* A is not a power of B*B the standard convergence theory for ill-posed 
problems cannot be used. Instead we apply the variable Hilbert scale theory and 
obtain convergence results for Eddington correction and Lorentz deconvolution of 
Gaussian and Voigt spectra in particular. Knowing these error bounds provides 
some insight into the choice of the enhancement operators B which goes beyond 
the range condition range(A) C range(B). 

By a change of perspective one interprets resolution enhancement as an appli- 
cation of an unbounded operator R. In the case of this paper, R = B^ 1 for the 
integral operator B. Another larger class of such enhancements is obtained when 
R is a differential operator. The theory of the application of such operators is 
covered in the recent book [ST] by Groetsch. A specific algorithm for numerical 
differentiation based on averaging and differences which converges with the size of 
the sampling with is analysed in [3J. The important question of the choice of the 
amount of differentiation for enhancement is discussed in [4] . 

If the broadening operator is known explicitly and is a convolution a different 
approach to resolution enhancement is based on the dilation (or rather contraction) 
of the spectral lines. Error bounds can also be obtained and a variant of variable 
Hilbert scales, the dilational Hilbert scales has been introduced to perform this 
analysis in [22 j . The approach has a particular appeal in practice as it does not 
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Figure 2. Regularised Lorentzian correction of a Gaussian 



introduce any satellite maxima. Such maxima might still occur, however, when 
data errors are large and regularisation has to be used. 

There is a substantial practical literature on separating overlapping line-shapes 
which cannot be covered here in any detail. As an example of a method which 
uses extra information, i.e., the ratio of the heights of two lines and the distance 
between them is the Rachinger correction formula [40]. This formula allows the 
determination of the corresponding line strengths Uj even without knowledge of the 
shapes a(-,Xi). In a sense, this is also what spectral enhancement methods attempt 
to achieve - but without any extra information. 

Related to the problem of spectral enhancement is the statistical problem of 
deconvolution of a density. Convergence rates have been found for several such 
problems in . These problems are often severely ill-posed and very slow conver- 
gence rates are obtained. The reason for this is that one can only assume that the 
underlying density is k times differentiable. While the authors did not use spectral 
theory nor the variable Hilbert scale interpolation inequality for their results one 
can obtain similar results with these more modern tools. This work has been contin- 
ued and practical estimators are discussed (also for less severely ill-posed problems) 
in [T5]. An interesting adaptive approach to these statistical problems is discussed 
in |46j where similar convergence results are obtained as in our discussion but using 
different techniques for analysis and different algorithms, see also [35, 28 . It would 
certainly be of interest to investigate these approaches from an ill-posed problem 
perspective using variable Hilbert scales. 

Maybe the most important limitation of the above discussion relates to the fact 
that all the operators occurring are convolutions. As outlined in the discussion 
of the models on broadening, the Doppler broadening is not a convolution and 
one can see that the operator may be factorised into diagonal operators and a 
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convolution. The next natural step would be to utilise norm equivalences (possibly 
using wavelets) with the variable Hilbert scale interpolation theory to deal with 
such more general source conditions. 
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